Prom Finite Nuclei to the Nuclear Liquid Drop: Leptodermous Expansion Based on 

the Self-consistent Mean-Field Theory 
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The parameters of the nuclear liquid drop model, such as the volume, surface, symmetry, and cur- 
vature constants, as well as bulk radii, are extracted from the non-relativistic and relativistic energy 
density functionals used in microscopic calculations for finite nuclei. The microscopic liquid drop 
energy, obtained self-consistently for a large sample of finite, spherical nuclei, has been expanded in 
terms of powers of A~ 1/:i (or inverse nuclear radius) and the isospin excess (or neutron-to-proton 
asymmetry). In order to perform a reliable extrapolation in the inverse radius, the calculations 
have been carried out for nuclei with huge numbers of nucleons, of the order of 10 6 . The Coulomb 
interaction has been ignored to be able to approach nuclei of arbitrary sizes and to avoid radial 
instabilities characteristic of systems with very large atomic numbers. The main contribution to 
the fluctuating part of the binding energy has been removed using the Green's function method to 
calculate the shell correction. The limitations of applying the leptodermous expansion to actual 
nuclei are discussed. While the leading terms in the macroscopic energy expansion can be extracted 
very precisely, the higher-order, isospin-dependent terms are prone to large uncertainties due to 
finite-size effects. 

PACS numbers: 21.30.Fe, 21.60.Jz, 24.10.Jv, 
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I. INTRODUCTION 

Bulk properties of atomic nuclei that depend in a 
smooth way on the numbers of nucleons have been tradi- 
tionally described in terms of macroscopic models, e.g., 
the liquid drop model (LDM) or the droplet model; for 
reviews, see Refs. 0,0,0, 0,0- These phenomenological 
models, often augmented by a shell correction which is 
calculated using average single-particle potentials, have 
been tuned up to describe nuclear bulk properties to a 
high precision. At the same time, these models can be in- 
terpreted in the language of the leptodermous expansion 
that sorts the various contributions to the binding en- 
ergy of finite nuclei in terms that have transparent phys- 
ical meaning, e.g., volume, surface, symmetry, curvature, 
and Coulomb energy. 

On the microscopic side, self-consistent mean-field 
models employing density-dependent effective interac- 
tions or energy-density functionals are nowadays com- 
monly used in nuclear structure modeling. The most 
prominent of these are the Skyrme-Hartree-Fock (SHF) 
method, the relativistic mean-field (RMF) approach (as 
well as their Bogoliubov extensions), and the Hartree- 
Fock-Bogoliubov method with the finite-range Gogny 



force; for a recent review, see Ref. 7}. These models 
rely on effective energy-density functionals with typically 
6-10 parameters adjusted phenomenologically that pro- 
vide a global description of all nuclei throughout the nu- 
clear chart (with exception, perhaps, of the lightest ones). 
While the parameters of these models can be organized 
and interpreted using the low-energy effective field the- 
ory of quantum chromo dynamics 0, 0, EJ , they do not 
have an immediate interpretation in terms of the total 
number of nucleons (or nuclear radius) and neutron ex- 
cess. For that reason, it is both convenient and practi- 
cal 1 1 1| to characterize nuclear energy-density functionals 
in terms of certain macroscopic parameters. The usual 
starting point is the limit of the homogenous infinite nu- 
clear matter, which is simple to compute and which de- 
fines the leading LDM characteristics, e.g., the volume 
energy, symmetry energy, or incompressibility. 

Nuclei have a pronounced surface; hence, a proper 
characterization of surface properties is crucial. This 
task, however, is far from easy. The usual means of 
characterizing surface properties of the energy functional 
is through semi-infinite nuclear matter having a planar 
surface zone (see Ref. an d references quoted therein). 
Early attempts used semi-classical approaches to circum- 
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vent the enormous complexity of self-consistent calcula- 
tions for the corresponding mean field (see, e.g., Ref. 
for the case of SHF). The limited self-consistent calcu- 
lations in SHF [3 HU and RMF [l6| indicate that the 
quantum effects are non-negligible to the extent that they 
change the surface parameters by about 5 %. 

In Ref. the volume and surface contributions to the 
energy density were extracted by assuming the Thomas- 
Fermi relation between the local density and the kinetic 
energy density. The authors concluded that, in the nu- 
clear surface zone, the gradient terms (absent in the ho- 
mogenous nuclear matter) are as important in defining 
the energy relations as those depending on the local den- 
sity. That is, the nuclear surface cannot simply be re- 
garded as a layer of nuclear matter at low density. 

In this work, we extract the macroscopic LDM param- 
eters by using a large sample of finite, spherical nuclei, in- 
cluding huge systems having 10 5 -10 6 nucleons. Based on 
the self-consistent SHF and RMF results, we extract the 
macroscopic information from the large-scale trends by 
subtracting fluctuating shell corrections. The Coulomb 
force is switched off to allow computation of very large 
nuclei. (We thus concentrate on the strong component 
of the nuclear interaction.) This way of analysis, using 
finite nuclei rather than semi-infinite matter, is concep- 
tually closer to existing nuclei, and it allows the deter- 
mination of curvature and surface symmetry effects. In 
essence, the aim of this paper is two-fold. First, by con- 
sidering finite, although huge, nuclei, we investigate the 
convergence of the macroscopic expansion. Second, by 
taking several SHF and RMF energy functionals, we ex- 
plore the relation between surface parameters and the 
nuclear matter features of the underlying forces. 

The paper is organized as follows. Section ITT1 discusses 
the macroscopic energy expansions. The details of our 
SHF and RMF models and shell-correction calculations 
are given in Sec. lIIII The extraction of LDM parameters 
is described in Sec. II VI they are discussed in Sec. \V\ 
Finally, conclusions are drawn in Sec. IVII 



II. THE MACROSCOPIC ENERGY 

According to the Strutinsky Energy Theorem [Tsl figt 
I20L l2l| , the energy per nucleon can be decomposed into 
an average part (smoothly depending on the number 
of nucleons) and the shell-correction term that fluctu- 
ates with particle number reflecting the non-uniformities 
(bunchiness) of the single-particle level distribution: 



g ^(smooth) _j_ j^(shell) 

A ~ - 



(1) 



Macroscopic models, such as LDM, deal with the smooth 
part. Consequently, in the following we concentrate on 
the average part of the binding energy: 



The macroscopic energy can be parametrized in many 
ways. However, by far the most successful macroscopic 
mass expressions are those rooted in the liquid-drop 
model (LDM) and in the droplet model. They are re- 
spectively outlined in Sec. Ill Al and Sec. Ill Bl below. 



A. Liquid-drop model 

The LDM parameterizes the binding energy of the nu- 
cleus (Z, N) in equilibrium. Instead of proton and neu- 
tron numbers, it is convenient to express the LDM energy 
through the mass number and neutron excess: 



A = N + Z 



I = 



N - Z 



N - 



(3) 



The macroscopic binding energy per nucleon can be ex- 
panded as 

£(LDM) _ ^(smooth) fj^ j\ 

— Q-vol ^surf-^ ^ CL CUTV A / 

+ a J 2 4- a J 2 4 -1 / 3 ( 4 ) 

1 Lt sym ± 1 u, ssy m ± 

+ a {2) I 4 

All the terms in Eq. (J2J have an immediate physical in- 
terpretation. The bulk energy is given by the volume 
energy a vo i, and changes with the neutron excess are ac- 
counted for by the symmetry-energy term a s 



^(smooth) _ ^^(shell) 

~ A 



(2) 



J 2 and 

by the second-order symmetry-energy term aiylnl 4 . The 
most important finite-size correction is the surface energy 
«surf^4 _1 ^ 3 , followed by more subtle trends in terms of the 
curvature energy a curv A~ 2 / 3 and the surface-symmetry 
energy a ssym I 2 A^ 1 / 3 . 

The sorting in columns indicates the level of impor- 
tance of the terms. Two different sorting criteria are used 
simultaneously: an expansion of finite size effects (=sur- 
face effects) in terms of powers of A~ x l 3 (proportional to 
inverse radii) and, parallel to it, an expansion in terms of 
the neutron-to-proton asymmetry I 2 . The second-order 
symmetry energy term cx I 4 is not always included in the 
macroscopic binding energy expression. It has been con- 
sidered, e.g., in the context of the Thomas- Fermi model 
[22| and in a discussion of strongly asymmetric matter 
within the RMF [2^ . We find that such a term appears 
naturally in the hierarchy of Eq. and we shall demon- 
strate that it is naturally present in the microscopic LDM 
expression. 

At this point, it is worth noting that the shell energy 
per nucleon, <5£( shc11 ), scales with mass as A~ 2 / 3 |24| . 
i.e., it has the same dependence on the nuclear radius 
as the curvature term. Consequently, uncertainties asso- 
ciated with the extraction of shell corrections from self- 
consistent results, and the presence of higher-order fluc- 
tuating terms that are not accounted for by the Strutin- 
sky procedure, can seriously impact the values of higher- 
order terms in the leptodermous expansion. We will see 
it very clearly in the results presented in Sec. IIVI 
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The ansatz |J2J does not include explicit information 
about the nuclear radius. In fact, the LDM tacitly as- 
sumes a fixed radius 



Ro = rsA 1 / 3 



(5) 



where the Wigner-Seitz radius is typically r s =1.14-1.20 
fm, which defines the saturation (equilibrium) density 



Pa 



(6) 



A more general ansatz which allows a determination of 
the radius is provided by the droplet model presented 
below. 



d and e, the relative deviation in the bulk of the density 
p from its nuclear matter value p$: 



R — Rq 
E = — R~ = 



P~ Po 
3p 



(7a) 
(7b) 



where Rq is the equilibrium radius of the droplet, see 
eq. J5j. Energy changes with d and e are considered 
explicitly. For instance, the volume term is augmented 
by a compression effect, 



ftvol 



Ovol 



(8) 



B. Droplet model 

The droplet model |2^| (see, e.g., Ref. Q for a recent 
implementation) includes the effect of the neutron skin 
and non- uniformities in the nuclear density. The two cru- 
cial parameters of the droplet model that describe devia- 
tions from the equilibrium arc the neutron skin thickness 



where the nuclear incompressibility coefficient is 



d 2 E 



K ^JP~-A 



(9) 



P=P0 



The droplet model binding energy per particle is con- 
sidered a function of A, I, e, and d: 



£(drop) = ^(drop)^^^ 

= avoi + cwA- 1 / 3 + d CUIV A- 2 / 3 + 2a surf A- 1 / 3 e + fe 2 (10) 
+ a sym I 2 + d ssym A- 1 / 3 f(I,d)-3a' sym I 2 e + di y ll 4 



where a' sym is defined as 



da s 



dp 



P=P0 



By substituting (|12fl in i|10|) , one arrives at the LDM ex- 
pression (@J where the leading parameters a vo \, a sur f, and 
(11) a sym remain unchanged while the higher-order parame- 
ters are redefined as: 



The function /(/, d) is assumed to be quadratic in / and 
d; it is determined by minimizing the energy with respect 
to <i, see below. 



C. Relation between droplet model and LDM 
expansions 

In this section, we recall the relation between the LDM 
expression (0J and the more detailed droplet mass for- 
mula (|10fl . The form (0} is valid at the energy in equi- 
librium (i.e., where radius and/or density have been ad- 
justed to minimize the energy for a given nucleus) while 
equation l|10|) allows for separate tuning of e and d. The 
equilibrium energy £ ( Ted >(A, I) is obtained by minimizing 
£( dr °p)(A, j £ ^ ^ w £th respect to e and d. At the equilib- 
rium, d = bl, where &ssl.4 fm. For the radial expansion 
e, one obtains 



^a^jA- 1 / 3 + 3a' sym I 2 
K 



O-curv 


— G-curv 


O 2 Po 
za surf j£ ! 


&ssym 


— ^ssym — 


oa sur f a sym 


sym 


- 5« - 

sym 


V ) 2 ^ 

2 V "sym 7 



(13a) 
(13b) 
(13c) 



It is seen that while the leading parameters, a vo i, a sur f, 
and o sym , are defined unambiguously in both LDM and 
droplet model expressions, the higher-order terms differ. 
We shall determine the parameters of the leptodermous 
expansion from the calculated ground state binding en- 
ergies. The corresponding mean-field configurations are 
stable points; hence, the equilibrium philosophy of the 
LDM should apply. To deduce the droplet model param- 
eters a curv , cissym, and a sym , one should use the relations 
( 12 ) 03- 
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TABLE I: LDM coefficients (in MeV) for the self-consistent mean-field models applied in this work. The first block (columns 
2 and 3) shows, for reasons of completeness, bulk droplet-model parameters po and K. The second block (columns 4-6) shows 
the bulk parameters of the LDM as computed in the limit of homogenous nuclear matter. The third block (semi-bulk; column 
7) shows surface parameters from semi-infinite nuclear matter calculations, where available (Ref. |Tsjl for SHF and Ref. |2(| for 
RMF). The results for the parameters deduced from finite nuclei as described below are shown in columns 8-10. The theoretical 
uncertainty on surface energies is 0.05 MeV for SHF and 0.1 MeV for RMF. The curvature energies are reliable within 0.5 
MeV. The surface-symmetry energies have an uncertainty of about 2 MeV for SHF and 10 MeV for RMF. The lowest rows 
show the LDM coefficients adjusted to data on finite nuclei. The corresponding parameters of the droplet model (1101 can be 
deduced from the given LDM parameters through the relations (I13H plus the empirically supported assumption H6aH . 
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semi- bulk 




from finite nuclei 




Model 


Po 


K 








(NM) 
a surf 


^surf 


ftcurv 


Issym 


SkM* 

SkP 

BSkl 

BSk6 

bLy4 

SLy6 

SkI3 

SkI4 

SkO 


0.1603 
0.1625 
0.1572 
0.1575 

n i c:o£ 

0.1590 
0.1577 
0.1601 
0.1605 


216.6 
201.0 
231.4 
229.2 
zoU.l 

230.0 
258.1 
247.9 
223.5 


-15.752 
-15.930 
-15.804 
-15.748 
-lo.y 1 1 
-15.920 
-15.962 
-15.925 
-15.835 


30.04 
30.01 
27.81 
28.00 

31.96 
34.84 
29.51 
31.98 


95.25 

40.43 

15.76 

35.67 

yo.y i 

99.48 

212.47 

125.80 

163.50 


17.70 
18.22 
17.54 

17.74 


17.6 
18.2 
17.5 
17.3 

lo.4 

17.7 
18.0 
17.7 
17.3 


9 

9.5 
9.5 
10 

n 

9 

10 

9 
9 
9 


-52 
-45 
-36 
-33 

-04 

-51 
-75 
-34 
-58 


NLl 
NL3 
NL-Z 
NL-Z2 


0.1518 
0.1482 
0.1509 
0.1510 


211.3 
271.7 
173.0 
172.4 


-16.425 
-16.242 
-16.187 
-16.067 


43.48 
37.40 
41.74 
39.03 


311.18 
269.16 
299.51 
281.40 


18.5 
17.7 


18.8 
18.6 
17.8 
17.4 


9 
7 
9 
10 


-110 

-86 
<-125 
-90 


LDM [4] 
LDM [5j 
LSD [5J 


0.153 
0.1417 
0.1324 




-16.00 
-15.848 
-15.492 


30.56 
29.28 
28.82 






21.1 
19.4 
17.0 


3.9 


-48.6 
-38.4 
-38.9 



III. SELF-CONSISTENT MODELS 

We employ two variants of self-consistent mean-field 
models: SHF and RMF. They are explained in detail in 
Ref. 0- Both approaches provide a functional form for 
the energy density with a good handful of free parame- 
ters. These have been adjusted to phenomenological data 
by different groups and with different bias. Thus there 
exist various parametrizations on the market which pro- 
vide a fairly good description of basic nuclear bulk prop- 
erties in the valley of stability, but differ in other aspects 
as, e.g., excitations, fission barriers, neutron matter prop- 
erties, or electromagnetic form factors. 

In this work, we have chosen a small subset of Skyrme 
forces which perform well for the basic ground-state prop- 
erties and have sufficiently different properties which al- 
lows one to explore the possible variations among param- 
eterizations. This subset contains: SkM* [27|, SkP [28| . 
SLy4 IH, SLy6 H3, SkO HU, BSkl H2, BSk6 [H, and 
Skll, SkI3, and SkM from Ref. 0. SkP, SkO, and BSkl 
have effective nucleon mass around one, leading to a com- 
paratively large density of single-particle levels. All other 



SHF forces employed here have smaller effective masses. 
Interesting here is the double BSkl with BSk6. Both 
forces were fitted using a similar strategy and data pool, 
but have different effective mass m* jm — 1.05 for BSkl 
and 0.8 for BSk6. The force SLy6 was adjusted with 
particular emphasis on isotopic trends and neutron mat- 
ter. The functionals SkI3, SkM, and SkO have a general- 
ized isovector spin-orbit interaction compared to all other 
forces. Some forces, i.e. SkP, BSkl, and BSk6, include a 
J 2 term (where J is the spin-orbit tensor) with a cou- 
pling constant related to the surface and effective mass 
terms, while all others do not. All the selected forces 
perform reasonably well concerning the total energy and 
radii for nuclei close to the valley of stability, with some 
different bias on particular observables. In particular, 
BSkl and BSk6 are fits to all available masses, but only 
masses. 

As seen in Table the values for the volume energy 
coefficient a vo i, saturation density po, incompressibity K, 
and surface energy coefficient a sur f are quite similar, with 
slight systematic differences between SHF and RMF that 
already have been noticed earlier; see Q and references 
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given therein. Obvious large variations occur for proper- 
ties which are not fixed precisely by nuclear matter and 
ground-state characteristics. 

As in SHF, there exist many RMF parameterizations 
which differ in details. For the purpose of the present 
study we choose the most successful (or most commonly 
used) ones: NL1 j^, NL-Z 0, NL-Z2 and NL3 
[38J. The parameterization NL1 is a fit of the RMF along 
the strategy of Ref. j^. The NL-Z parametrization is 
a refit of NL1 where the correction for spurious center- 
of-mass motion is calculated from the actual many-body 
wave function, while NL-Z2 is a recent variant of NL-Z 
with an improved isospin dependence. The force NL3 
stems from a fit including exotic nuclei, neutron radii, 
and information on giant resonances. All the above pa- 
rameterizations provide a good description of binding 
energies, charge radii, and surface thicknesses of stable 
spherical nuclei with the same overall quality as the SHF 
model. As seen in Table [I] however, the nuclear matter 
properties of the RMF forces show some systematic dif- 
ferences as compared to Skyrme forces. All RMF forces 
have comparable small effective masses around m* /m w 
0.6. (Note that the effective mass in RMF depends on 
momentum; hence the effective mass at the Fermi energy 
is approximately 10% larger.) Compared with SHF mod- 
els, the absolute value of the energy per nucleon is sys- 
tematically larger, with values around —16.3 MeV, while 
the saturation density is always slightly smaller with typ- 
ical values around 0.15 nucleons/fm 3 . The incompress- 
ibility of the RMF forces ranges from low values around 
170 MeV for NL-Z to K=270 MeV for NL3. There are 
also differences in isovector properties; the symmetry en- 
ergy coefficient of all RMF forces is systematically larger 
than for SHF interactions, with values between 37.4 MeV 
for NL3 and 43.5 MeV for NL1. 

The absolute variation of the LDM energy coefficients 
found in Table [I] should be put into perspective. There 
is a notable difference between the variation of the coef- 
ficients in the LDM expansion on the one hand and the 
actual variation of the energy on the other hand. Let 
us consider the heavy nucleus 250 Fm with A = 250 and 
/ = 0.2 as an example. The difference in a vo i between 
SkM* and NL1 is 0.67 MeV, which seems to be small. 
It leads, however, to an energy difference of about 160 
MeV, which amounts to a significant fraction (~10 %) of 
the total binding energy of this nucleus. The difference in 
the symmetry energy coefficient, a sym between the same 
two forces is 3.44 MeV and appears to be much more 
significant that the difference in a vo \. However, the I 2 
factor suppresses the difference in the symmetry energy 
between SkM* and NL1 to 34.4 MeV. Even the 60 MeV 
difference between the surface symmetry energy coeffi- 
cients of SkM* and NL1 only gives about 100 MeV dif- 
ference in binding energy between the two forces, which 
is comparable with, but still smaller than, the energy 
difference arising from the small difference in the volume 
energy coefficient. 

It is interesting to check to what extent the differences 



in po, K, and a' syrn influence the macroscopic parameters 
and whether a simple correlation between spectroscopy 
and macroscopy can be found. This will be done in 
Sec. lVBl where we will consider a larger variety of param- 
eterizations and perform dedicated variations of special 
features as, e.g., the effective nucleon mass. 

According to Eq. J5J), the average energy is obtained 
by subtracting the shell-correction energy from the self- 
consistent value. The shell energy is computed using the 
same prescription as outlined in Rcfs. [4(J,[4l|]. This pro- 
cedure is based on the Green's function approach to the 
level density and the generalized plateau condition of Ref. 
|4^ . The advantage of this procedure is that it prop- 
erly takes care of the continuum positive-energy states 
which unavoidably come into play in the self-consistent 
approach. In our calculations, we include a large space 
of single-particle states up to 60 MeV above the Fermi 
energy. Since most of these states are continuum states, 
the contribution from a particle gas (treated in the same 
numerical box) has to be removed. To meet the general- 
ized plateau condition, it is assumed that the deviation 
of the smoothed level density from linearity is minimal 
in a wide energy interval [-50,-20] MeV. 

Our calculations are restricted to spherical symme- 
try; they were carried out using numerical techniques 
described in detail in Ref. 43]. The Coulomb force 
is ignored to allow an extension to nuclei of arbitrary 
sizes. Pairing correlations are neglected as well. How- 
ever, open-shell nuclei were treated in the filling approxi- 
mation in which we use a very small pairing force (factor 
of 10 smaller than usual). The center-of-mass (cm.) cor- 
rection is included. We take care to use precisely the cm. 
recipe that is attached to a given force • This is crucial 
because it is known that the actual form of the cm. cor- 
rection has a significant impact on the surface properties 



IV. EXTRACTION OF LDM PARAMETERS 
A. Bulk parameters 

The bulk parameters in the leptodermous expansion 
are those proportional to ^4°. They represent terms which 
do not vanish in the A — > oo limit. In the LDM ansatz 
(|4*|) the bulk parameters are: the volume energy constant 
a vo i, the symmetry energy constant a sym , and the second- 

(2) 

order symmetry energy parameter asym- The droplet 
model (|10|l additionally introduces the incompressibil- 
ity K, the density-slope of the symmetry energy a' sym , 
and the equilibrium density pq. All these parameters can 
easily be computed in the limit of the homogenous bulk 
nuclear matter, see, e.g., Ref. Q. 

Asymmetric nuclear matter shows an interesting trend, 
which sheds some light on the relation between the LDM 
and droplet model. The bulk part of the LDM energy 
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FIG. 1: Solid line: the effective symmetry energy parameter 
l|15fl in asymmetric nuclear matter computed with NL-Z2 as 
a function of neutron excess /. Dotted line: the trend d!6bl> 
using the bulk parameters of NL-Z2, see Table U] 



FIG. 2: The binding energy per nucleon £ (A, I) for isospin- 
symmetric (I = 0) nuclei computed with SLy6 as a function 
A _1,/3 before (solid line) and after (dotted line) subtracting 
the shell correction ^f' 8 ' 1011 - 1 . The fine dotted line connects the 
last value from finite nuclei with the nuclear matter limit. 



can be written as: 



£(LDM) = flvoi + flsym/ 2 + a (2) m/ 4^ (w) 



In order to concentrate on the isospin dependence, we 
introduce an effective symmetry energy parameter as 



This is consistent with Eq. I|16b|) and droplet-model pa- 
rameters displayed in Table [i] Finally, corroborating ev- 
idence of a very small ct sym term comes from Brueckner- 
Hartree-Fock studies of asymmetric matter based on re- 
alistic nucleon- nucleon and three-body forces [i^lliri . 



i off 

sym 



£(I)-£(I = 0) 
I 2 



(15) 



+ a (2) I 2 

1 sym 



The first line serves as a general definition. The second 
line then is specific to the bulk limit. 

Figure ^ shows the effective symmetry energy parame- 
ter obtained in the NL-Z2 model. Let us now recall that 

(2) (2) 

dsym is related to the droplet parameter as ym through 
Eq. I|13cp . This suggests that there is no explicit second- 
order isospin correction to the symmetry energy in the 
droplet model: 



",(2) 



0, 



which yields 



2 off 

sym 



(A=oo, J) 



-fa' ) 2 ^I 



A sym 



(16a) 



(16b) 



It is seen in Fig. ^ that the estimate <|16b|l matches the 
exact result extremely well. The deviations are small and 
predominantly cx I 4 , thus going beyond the present ex- 
pansion. (Higher-order isospin corrections were in fact 
considered in Ref. [2j| in the context of RMF and infi- 
nite nuclear matter.) We have checked that the assump- 
tion is well fulfilled for all the RMF and SHF func- 
tionals used in this work. The isospin dependence of 
a°y m (A = oo, /) in SHF is much weaker than in RMF. 



B. Isospin-independent surface parameters 

The leading isospin-independent parameters character- 
izing finite-size (surface) terms in the LDM are the sur- 
face and curvature energy coefficients. We deduce them 
from the systematics of binding energies of spherical nu- 
clei. 

Figure |21 shows the systematics of binding energies per 
nucleon predicted by SLy6. The smooth component ob- 
tained by subtracting <5£( shc11 ) is also indicated. At very 
large values of A (i.e. small A^ 1 / 3 ), the binding en- 
ergy per nucleon nicely converges to a straight line which 
demonstrates the validity of the hierarchical LDM (or 
droplet) ansatz and, not surprisingly, hints at the domi- 
nating role played by the surface energy as first leading 
correction to the bulk. The plot also illustrates the effect 
of quantum shell fluctuations. The actual binding en- 
ergy oscillates around the average trend; the amplitude 
of shell oscillations increases in lighter nuclei, consistent 
with the expected A _2 / 3 -dependence discussed earlier in 
Sec. lII Al By subtracting the shell correction, one obtains 
a fairly smooth trend, at least at this level of analysis. 

To extract the surface- and curvature-energy coeffi- 
cients, it is convenient to introduce the effective surface- 
energy coefficient: 



,(eff) 



,A 



-1/3 



(17a) 
(17b) 
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FIG. 3: Effective surface-energy coefficient If 7aH versus 
A~ 1/3 for several SHF (top) and RMF (bottom) function- 
al, as indicated. A straight-line fit to the data is marked by 
a dashed line. 

which is a function of system size A. The surface-energy 
coefficient <z sur f is obtained by extrapolating A^ 1 / 3, — > 0. 
The curvature a CU rv is then obtained from the slope of 
(tTTB . 

Figure[3]shows the effective surface energies. Note that 
I 



a sym(A-0 — 



Figure 01 demonstrates how such methodology works. 
The symmetry energy involves a difference of smoothed 
energies; hence, a difference of shell corrections. These 
corrections are prone to uncertainties, as we have already 
seen in Sec. lIVBl At small values of J 2 , the remaining un- 
controlled energy fluctuations are amplified in the finite 
difference (|18[1 . and the result is not reliable. Fortunately, 
at larger values of I > 0.1, where the isospin-dependent 
terms dominate over remaining shell fluctuations, one al- 
ways obtains a stable and well-defined plateau. The value 
of / is limited from above by the neutron drip line. In- 



these are drawn from smoothed energies (i.e. after sub- 
traction of the shell corrections) . The remaining fluctu- 
ations are due to higher-order shell effects [2£J, |2l| which 
cannot be accounted for by the generalized Strutinsky 
procedure. The construction (|17a|l of the effective sur- 
face energy amplifies those residual fluctuations dramat- 
ically; it is only by virtue of the smoothed energy that 
one can see any clear trend. Thus atop these remaining 
shell fluctuations, one can recognize a definite slope (note 
the very narrow energy window), which is related to the 
curvature energy. By performing a straight-line fit to the 
data, surface- and curvature-energy coefficients can be 
extracted. 

Note that different A^ 1 / 3, scales are used for SHF and 
RMF. The reason is that in SHF cases we have been able 
to carry up calculations for really huge nuclei (A = 10 6 ). 
This extends the scale to smaller A" 1 / 3 and allows us to 
ignore some data points for lighter systems where resid- 
ual shell fluctuations are large. Unfortunately, we were 
unable to approach similarly large nuclei in RMF; hence, 
we had to scale up to A" 1 / 3 = 0.2 to get sufficient data 
for extrapolation. In any case, one sees that one can 
extract reliably well surface and curvature energies from 
the trends displayed in Fig. |3| We estimate an uncer- 
tainty in a sur f to be about 0.05 MeV for SHF and 0.1 
MeV for RMF. The curvature coefficient is determined 
within about ±0.5 MeV. 



C. Surface-symmetry coefficient 

In order to deduce the isospin-dependent surface- 
symmetry coefficient, we come back to the effective sym- 
metry parameter 1)15( 1 now considering finite A and iden- 
tifying £ «-> £( smooth ). As it is convenient to subtract the 
known infinite-matter trend of Eq. (|16b|l . we introduce 
the reduced effective symmetry-energy coefficient: 



(18) 

I 

deed, at / w 0.3 the neutron Fermi energy becomes pos- 
itive and the self-consistent solution can no longer be 
trusted. Therefore, for further analysis, we introduce 
an /-averaged reduced effective symmetry-energy coef- 
ficient: 

= / dIalf m (A,I) I [ dl . (19) 
Jo.i I Jo.i 

The surface symmetry energy coefficient is obtained by 
plotting versus A~ 1 / z . The slope for small A -1 / 3 
corresponds to a ssym , similarly as it was done in Sec. lIVBl 



J_ ^(smooth) _ £ 



(smooth) i a n \ i 9( a sym) p , 



8 




FIG. 4: Top: reduced effective symmetry-energy coefficient 1181 1 in SLy6 versus I for two system sizes (t4=2840, left; A=80,000, 
right). Bottom: corresponding proton and neutron Fermi energies. 



for deducing the surface and curvature parameters. An 
effective surface-symmetry constant can thus be obtained 
from 



(eff) 
ssym 



( „rod _ „ ^ 4l/3 
I sym "sym I 



(20) 



Figure shows typical results of such analysis. It is 
gratifying to see that the effective symmetry-energy co- 
efficients (|19|) are consistent with the corresponding bulk 
values when extrapolating back to A -1 / 3 — > 0. The qual- 
ity of the deduced slope (— surface asymmetry) can be 
assessed by inspecting the rescaled quantity (|2*U)l shown 
in the upper panel. Clearly, very large A, i.e., very small 
A -1 / 3 , are required to see convergence to a strictly hori- 
zontal trend. With the present data set, one may attach 
~10 % relative uncertainty (or about 2 MeV absolute 
error) to a ssym . 



D. Radii 

By employing the equilibrium value of e <|12|) . we can 
estimate the droplet-model radius as 



R = R a (l-e). 



(21) 



It is worth checking the performance of that recipe. To 
this end, we have extracted diffraction radii and r.m.s. 
radii from the SLy6 calculations of large (and huge) 
spherical nuclei with N—Z. Figure [5] displays the nu- 
clear radii corrected for shell fluctuations up to A=4000. 
It is seen that the estimate (|21[l evaluated with the SLy6 
parameters of Table [fl nicely approximates the actual re- 
sults. The r.m.s. radii are systematically larger than the 
diffraction radii which is no surprise because in the Helm 
model both are related via 



= \~\R. 



Miff 



5a 2 



(22) 



where a is the surface thickness coefficient (related to the 
first maximum of the form factor) with a typical value of 
1 fm. The estimate <|IVD(1 is indicated in Fig. and it 
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FIG. 5: Bottom: the averaged reduced symmetry-energy 
coefficient l|19jl in SLy6 and NL-Z2; Top: the corresponding 
effective surface-symmetry coefficient 1201 , 



is shown to be fully consistent with the self-consistent re- 
sults for the r.m.s. radii. Moreover, our analysis makes it 
very clear that the diffraction radius (or box-equivalent) 
is indeed the appropriate quantity entering the droplet 
model. 



V. RESULTS AND DISCUSSION 
A. Existing energy functionals 

Table|Ucollects the calculated LDM parameters for the 
SHF and RMF parametrizations introduced in Sec. II I II 
The bulk parameters are computed for the homogenous 
nuclear matter. The surface-related parameters are de- 
duced from the present analysis of finite nuclei. We 
also show for comparison the surface-energy coefficients, 

(NM) 



^surf ' obtained from quantum mechanical calculations 
for the semi- infinite nuclear matter [TR I2r1 | . The com- 
parison of c^urf^ with the values of a sur f deduced from 
finite nuclei shows a nice agreement between both meth- 
ods of analysis. 

The surface-energy and curvature-energy coefficients 
are fairly robust quantities in that the variations through- 
out all forces are small. The curvature energy coefficient, 
9.5 ± 0.5 MeV, is smaller than other LDM constants but 




0.05 0.10 0.15 0.20 

.-1/3 



FIG. 6: Nuclear radii (multiplied by A~ 1/:i to remove the 
overall mass dependence) for N — Z nuclei computed with 
SLy6 and (partly) corrected for shell fluctuations using the 
Strutinsky method. The lower curve shows the diffraction 
radii and the upper curve shows the r.m.s. radii ^/5/3r rms . 
The nuclear matter Wigner-Seitz value of 1.145 mi is indi- 
cated by a box and a horizontal dashed line. The droplet 
model estimate 1211 of the diffraction radius and the esti- 
mate (11V Dl of the r.m.s. radius are shown (dotted lines). It 
is gratifying to see the agreement between the droplet-model 
estimates and the results of self-consistent calculations for fi- 
nite nuclei. 



certainly nonzero. Moreover, its variations remain small 
throughout all forces "from the shelves", below the un- 
certainties of our analysis. In that context, we want 
to comment on the investigations of Ref. [4£| (see also 
discussion in [j|) where it was claimed that one needs 
small curvature-energy coefficients to obtain a good fit 
to fission barriers. While we find about the same cur- 
vature energy for a variety of forces, the fission barriers 
of actinides and superheavy elements investigated in |49j 
turned out to be strongly force-dependent for a similar 
pool of interactions. We suspect that another key param- 
eter affecting fission barriers lies in the isovector channel. 
As one hint we will see in Sec. IV Bl that strong variations 
of the symmetry-energy coefficient can influence a C urv 

The surface energy is more important and its (small) 
variations are larger than the uncertainty of about ±0.1 
MeV. The observed trends can be sorted in various ways. 
One important aspect is that a sur f depends to some ex- 
tent on the strategy used when fitting a functional. It is 
obvious that the weight which was given to light nuclei in 
the fit has an influence on the surface coefficients, as well 
as whether surface properties were considered (e.g., the 
surface thickness in the case of SkI3 and SkI4). Another 
crucial ingredient in the trends with A" 1 / 3 is the way 
in which the center-of-mass (cm.) correction was im- 
plemented 0| . In the sample of forces considered here, 
there are two forces, SLy4 and SLy6, which were fitted 
with precisely the same strategy but differ in their treat- 
ment of cm. correction. The difference of 0.7 MeV in 
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flsurf is quite remarkable and represents mainly the dif- 
ference between the recipes for the cm. correction used 
We think that the recipe used in SLy6, namely to 
compute the cm. energy from the given mean field state, 
is better microscopically motivated. Thus the lower value 
for a S urf is probably more realistic. However, this has yet 
to be explored in more detail. 

The situation for all the isospin-dependent LDM pa- 
rameters is quite different. That begins with the large 
discrepancy between a sym -values for RMF and SHF. In 
fact, even within SHF alone there is a much larger varia- 
tion in the symmetry energy than appears from Table 
see, e.g., the discussion in Ref. 50]. Extended RMF 
functionals also show significant uncertainty in a sym [5l| . 
Even more pronounced are fluctuations in the surface- 
symmetry coefficient o ssym . By inspecting Table [i] one 
can see a rough correlation between a sym and a ssym . This 
will be put on a firmer ground in Sec. IV Bl below where 
systematic variations of functionals are discussed. 



B. Systematic variations of functionals 

The discussion of Table [I] in Sec. IV Al has indicated 
several features which deserve closer inspection. To this 
end, we perform systematic variations of key properties 
of the functional, as, e.g. effective nucleon mass to* /to or 
symmetry energy a sym . The strategy is to vary only one 
chosen property while keeping all other features fixed. A 
set of SHF forces has been produced that way by fitting 
the parameters always to the same set of data (energies, 
charge form factors, spin-orbit splitting) while putting a 
constraint on the required additional feature. This was 
done formerly for the purpose of studying trends in the 
giant resonances [5(| . We consider these families of func- 
tionals in this work to inspect trends and correlations in 
a systematic manner. 

Figure [7J shows results of two sets of calculations. The 
first group in the left panel contains functionals with sys- 
tematically varied effective mass to* jm. The functionals 
belonging to the second set, shown in the right panel, 
vary a sym . Looking at the left panel, we find sizable 
variations with effective mass concerning a sur f, K, and 
a' synl . All those parameters decrease with m*/m. The 
surface-symmetry coefficient a ssym is fairly insensitive to 
to* /to. The curvature energy coefficient a curv (not shown 
here) varies only between 9 and 10 MeV which is very 
little in view of the theoretical uncertainty in this pa- 
rameter. Different trends are seen when the variation 
in a sym is considered (Fig. [7J right). There is a dra- 
matic change of the isospin-dependent terms, namely the 
density-dependent symmetry-energy coefficient a' sym and 
the surface-symmetry coefficient a ssym . The magnitude 
of variations in a sur f is comparable to that from the first 
set. That is, a lower value of G S ym can compensate for a 
larger effective mass. (A similar conclusion, in a slightly 
different context, has been drawn in Ref. 52].) And here 
is the first time that we see a handle on the curvature 



coefficient a curv : it weakly decreases with a sym . 

The above trends explain most of the results for the 
standard Skyrme functionals displayed in Table [I] How- 
ever, the comparison between BSkl and BSk6 leaves 
some puzzles. The functional BSk6 has a lower effective 
mass than BSkl (0.8 versus 1.05), and yet, the surface- 
energy coefficient slightly shrinks. This can probably be 
related to an improved treatment of the center-of-mass 
correction in BSk6. 

As already mentioned, there exist different prescrip- 
tions for defining the spin-orbit interaction in SHF Q. 
One variation concerns the possible contribution from 
the kinetic forces resulting in a J 2 -term in the energy 
functional (J is the spin-orbit tensor). We have studied 
this point within our systematic calculations. The results 
in Fig. [7] employ a version without (left panel) and with 
(right) that term. The same conditions are met for the 
left set at to*/to = 0.81 and the right set at a sym = 32 
MeV. All quantities shown are insensitive to that change, 
except for the density-dependent symmetry energy coef- 
ficient o.g ym , which exhibits a surprisingly large sensitiv- 
ity to J 2 . We have checked this point in more detail 
and concluded that the spin-orbit tensor term basically 
changes the offset of a' sym but has no influence on its 
trends. Anyway, it is noteworthy that a change of shell 
structure (here, via the spin-orbit interaction) can have 
such a dramatic effect on a bulk property of the func- 
tional. 



VI. CONCLUSIONS 

In this study, we present a systematic survey of nuclear 
surface properties in terms of the liquid drop model. Sur- 
face, surface-symmetry, and curvature energy coefficients 
are deduced as they are defined in the LDM, namely from 
the trends oc A -1 / 3 and oc A~ 2 I 3 in the binding energies 
of finite nuclei over a wide range of sizes. In order to 
achieve sufficiently precise values, we have evaluated a 
smooth background of binding energies by subtracting 
the shell corrections, and we considered huge nuclei con- 
taining up to 10 5 nucleons. 

Our calculations show how the bulk-matter limit is re- 
covered in finite nuclei. While it has been known from 
earlier studies [5^] that semi-classical features arc re- 
vealed only for nuclei with A > 5000, we found that ex- 
tremely massive nuclei are essential in order to pin down 
unambiguously the macroscopic surface-related parame- 
ters. The question emerges what role the LDM back- 
ground plays for actual nuclei which are extremely small 
at that scale. 

What is the influence of uncontrolled residual shell ef- 
fects on LDM parameters when only dealing with small 
sample of actual nuclei? The recent SHF work [54| that 
used a sample of "small" nuclei to extract the symmetry- 
energy and surface-symmetry energy coefficients can pro- 
vide a hint. While in some cases their results for a ssvm 
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FIG. 7: Dependence of LDM parameters on effective mass (left) and the symmetry-energy coefficient (right). The incompress- 
ibility is not shown in the right panel because it does not change. The curvature energy coefficient is not shown in the left 
panel because it varies very little as a function of m*/m. The extended Skyrme functionals of Ref. |H(i|| are used. 



are close to ours (e.g., they obtain a ssym =-49.2 MeV for 
SkM* versus -52 MeV here), there are forces for which 
the difference is fairly large (e.g., they obtain a ssym =- 
49 MeV for SkO while we get -58 MeV), and no clear 
tendency can be observed. 

In fact, the importance of shell effects depends on the 
observable. Energies as such exhibit a quick quenching of 
shell effects with increasing A (see Fig. EI) - However, en- 
ergy differences are required when evaluating the effective 
surface energy. As illustrated in Fig. |3 these differences 
reveal shell effects up to any size. This annoying fea- 
ture limits the precision with which one can deduce the 
higher-order LDM parameters, in particular the isovec- 
tor parameters which are roughly determined empirically. 
Those are also much harder to pin down in the present 
analysis, and probably in any analysis, as the available 
span of I in bound nuclear systems is rather small. 

The leading LDM parameters can be determined suf- 
ficiently well as to allow a thorough analysis. First, 
we have studied a broad selection of widely used "stan- 
dard" SHF and RMF parameterizations. We find that all 
isoscalar energy coefficients, including surface and cur- 



vature ones, show only a few percent variation, while 
the isovector energy coefficients might differ by a fac- 
tor of two or even more between forces. In the second 
step, we have worked out some interrelations by a sys- 
tematic variation of forces. We reconfirm that isovector 
features are much more sensitive to parameter changes 
than the isoscalar terms. We have also spotted a sur- 
prisingly strong interrelation between the spin-orbit force 
and some LDM properties (the density dependent asym- 
metry energy coefficient). This demonstrates that the 
shell structure and smooth LDM background are inti- 
mately connected. The variation of the LDM coefficients 
does not, however, directly translate into similar varia- 
tions of the corresponding energies. 

Finally, we note that in this work we follow a strictly 
"empirical" approach relating (shcll-corrcctcd) binding 
energies to the LDM parameters. One can amplify (and 
extract) some of the LDM parameters better by using 
other, more sensitive, observables, such as fission barriers 
or energies of superdeformed states. Such studies are 
currently being pursued. 
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